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We investigate the transport properties in a zigzag silicene nanoribbon in the presence of an external electric 
field. The staggered sublattice potential and two kinds of Rashba spin-orbit couplings can be induced by the 
external electric field due to the buckled structure of the silicene. A bulk gap is opened by the staggered 
potential and gapless edge states appear in the gap by tuning the two kinds of Rashba spin-orbit couplings 
properly. Furthermore, the gapless edge states are spin-filtered and are insensitive to the non-magnetic 
disorder. These results prove that the quantum spin Hall effect can be induced by an external electric field 
in silicene, which may have certain practical significance in applications for future spintronics device. 
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Recently, the quantum spin Hall effect (QSHE) has 
attracted significant interests in the fields of condensed 
matter physics and material science as it constitutes a 
new phase of matter and has potential applications in 
spintronics.^"^ The novel electronic state with a nontriv- 
ial topological property and time-reversal invariance has 
a bulk energy gap separating the valence and conduc- 
tion bands and a pair of gapless spin-filtered edge states 
sat the sample boundaries. The QSHE has been first 
predicted by Kane and Mele in graphene in which the in- 
trinsic spin-orbit coupling opens a band gap at the Dirac 
points.^ However, the QSHE can occur in graphene only 
at unrealistically low temperatures since the intrinsic 
spin-orbit coupling in graphene is rather weak.^"^ There- 
fore, it is crucial to search new materials with strong 
spin-orbit coupling for realizing the QSHE. Recent the- 
ories and experiments provide evidences of the QSHE in 
two-dimensional HgTe-CdTe quantum wells. ^''^ 

Very recently, a close relative of graphene, a slightly 
buckled honeycomb lattice of Si atoms called silicene 
has been synthesized. ^"^'^ Silicene can be well compatible 
with current silicon based electronic technology. Many 
progresses in the study of silicene have been made, both 
experimentally and theoretically. For example, electronic 
properties and the giant magnetoresistance in silicene 
have been reported. ^^"^^ Moreover, almost every striking 
property of graphene could be transferred to silicene. ^'^'^^ 
It has been theoretically shown that the strong intrin- 
sic spin-orbit coupling in silicene may lead to detectable 
QSHE.5'19-23 

In this paper, we provide systematic investigations on 
the band structures and electron transport properties of 



silicene in the presence of an external electric field. Sil- 
icene consists of a buckled honeycomb lattice of silicon 
atoms with two sublattices A and B. We take a silicene 
sheet on the x ^ y plane, and apply the electric field in 
z direction. The electric field generates a staggered sub- 
lattice potential between silicon atoms at A sites and B 
sites due to the buckled structure of the silicene. On the 
other hand, two kinds of Rashba spin-orbit coupling, re- 
ferring to the nearest and next-nearest neighbor hoppings 
respectively, can also be tuned by the external electric 
field. We find that a gap can be opened by the staggered 
sublattice potential and gapless edge states are induced 
in the gap by Rashba spin-orbit coupling. We predict 
that the QSHE can be observed by applying an exter- 
nal electric field in silicene even if the intrinsic spin-orbit 
coupling in the system is very weak. 

In the tight-binding representation, the silicene sample 
with an external electric field can be described by the the 
following Hamiltonian:^° 

{ij)a iOL 



(1) 



where creates an electron with spin polarization a at 
site i; (ij) and {{ij)) run over all the nearest and next- 
nearest neighbor hopping sites, respectively. The first 
term is the nearest-neighbor hopping with the transfer 
energy t = 1.6eV. The second term is the staggered 
sublattice potential term, where /Xy = ±1 for the A 
(B) site and Si is the potential energy induced by the 
external electric field. The third and fourth terms, re- 
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spectively, represent the first Rashba spin-orbit coupling 
associated with the nearest neighbor hopping and the sec- 
ond Rashba spin-orbit couphng associated with the next- 
nearest neighbor hopping. Both of them are induced by 
the external electric field. Here <j = {<7x,<yy,<yz) is the 
Pauli matrix of spin and cH^j = dij/\ d^j \ with the vector 

dij connecting two sites i and j. The intrinsic SOC term 
has been ignored intentionally since the main focus of 
this work is the Rashba terms and the staggered poten- 
tial, which can be tuned by the external electric field. 

We assume that the temperature is set to zero and 
two semi-infinite silicene ribbons are employed as left and 
right leads. The two-terminal conductance of the system 
can be calculated by the nonequilibrium Green's function 
method and Landauer-Biittiker formula as 

G{E) = ^TT[TL{E)G^E)rHiE)G^{E)], (2) 

where Tp{E) = i['El,{E) - i:^{E)] is the line-width func- 
tion and G''{E) = [G''{E)]'< 1/[E - U^en - S£ - 
is the retarded Green function with the Hamiltonian in 
the center region Hcen-^"' The self-energy due to the 
semi- infinite lead-p can be calculated numerically.^^ 

With the help of the nonequilibrium Green's function 
method, the local current flowing on site i with spin a 
can be expressed as 

ja' ja' 
- Gf„,ja'{^^'t)Hja' Aa\^ (3) 

where Gf^ = i{Cj^,Cia-) is the matrix element of the 
lesser Greens function of the scattering region and Jia,ja' 
is the current from site i to j. After taking Fourier trans- 
formation, the local current Jia-ja' can be calculated as 

J,.,.' = / ^MHr.,.'G<,^.,^m- (4) 

Eq. (4) has been widely used in the local-current studies 
of tight-binding models. ^^"^^ When the sample is at zero 
temperature and the applied voltage is small, by applying 
the Keldysh equation G< = G''(iri/i + irij/ij)G" with 
the Fermi distribution function fp = foiE + eVp),'^^ the 
Eq. (4) can be written as 

Jraja' = ^ / dElm{H,„,,^, [G^ {T L + rfl)G"],,.,,4 

" J —CO 

+ Vn)lm[H,,^,,,G],,^,AE)], (5) 

where Vl and Vr are the voltages at the Lead-L and 
R, respectively. G"(£;) = G''iE)rL{E)G''iE) is electron 
correlation function. The first part of Eq. (5) can only 
generate the equilibrium current and does not contribute 
to the transport, so it can be dropped out in present 
work. It is the second part that gives rise to the nonequi- 
librium current. 
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FIG. 1. Calculated energy bands in the zigzag nanoribbon 
geometry for (a) Si = 0, A^i = and \r2 = 0, (b) Ei = O.lt, 
\ri = and Xr2 = 0, (c) £i = O.lt, Xm = l.Ot and Xr2 = 
0.2t, and (d) Si = O.lt, Xri = 0.05t and Xr2 = 0.6t. 



In the following numerical calculations, we use the hop- 
ping energy t as the energy unit. The width of the zigzag 
ribbon is 59a, where a is the silicon-silicon distance. In 
Fig. 1 we show the energy bands obtained from diago- 
nalizing the tight-binding Hamiltonian (1) with various 
parameters for a zigzag nanoribbon. The nanoribbon 
with only the nearest-neighbor hopping shows a semi- 
metallic behavior, as shown in Fig. 1 (a). An energy 
gap can be opened due to the inversion symmetry break- 
ing induced by the staggered sublattice potential and the 
magnitude of the gap is 2£i (see Fig. 1 (b)). When the 
first and second Rashba spin-orbit couplings induced by 
the external electric field are taken into account properly, 
which turn silicene from normal insulating to quantum 
spin Hall regime, gapless edge states appear within the 
band gap (see Figs. 1 (c) and (d)). The gapless edge 
states with different spins connect the conduction band 
and valence band. As usual, these gapless edge states are 
originated from the nontrivial topological orders in the 
bulk. According to different values of the first Rashba 
spin-orbit coupling, the QSHE induced by the external 
electric field can be divided into two types, QSHEl and 
QSHE2. In QSHEl, when the first Rashba spin-orbit 
coupling is strong, the edge states traverse the bulk gap 
within each valley, as shown in Fig. 1 (c). However, 
in QSHE2 when the first Rashba spin-orbit coupling is 
weak, the edge states inter-connect two valleys. More- 
over, they bend and give rise to "subgaps" around K and 
K' , which makes the structure of propagating channels 
complicated in the bulk gap. 

To investigate the QSHE induced by the external elec- 
tric field in more details, the configurations of the spin- 
dependent local-current-flow vector are plotted in Fig. 
2. We focus only on the left-injected current. The Fermi 
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FIG. 2. The local spin-up (al and bl) and spin-down (a2 
and b2) currents density vector distribution for (a) Si = O.li, 
Am = l.Ot and Xr2 = 0.2t (QSHEl) with E = 0.05t, and 
for (b) Ei = O.lt, \ri = 0.05t and \r2 = 0.6t (QSHE2) with 
E = 0.005t. 



vector in such regions are plotted in Fig. 3. We find that 
the edge states are not fully spin-filtered when they are 
inside the bulk band. However, in this case, the electrons 
flow along the lower edge (see Fig. 3 (a)), while the holes 
flow forward along the upper edge (see Fig. 3 (b)). 




energy is set to be £; = O.OSt for QSHEl (Fig. 2 (a)) 
and E = O.OOSt for QSHE2 (Fig. 2 (b)). For these Fermi 
energies, there arc only the lowest transmission channels, 
i.e., the gapless edge states. For QSHEl, the spin- up lo- 
cal currents locate mainly on the lower edge (see Fig. 2 
(al))and the spin-down local currents locate mainly on 
the upper edge (see Fig. 2 (a2)). Contrary to QSHEl, 
for QSHE2, the spin- up local currents locate mainly on 
the upper edge (see Fig. 2 (bl)) and the spin-down lo- 
cal currents locate mainly on the lower edge (see Fig. 2 
(b2)). These results show that the gapless edge states 
are spin-filtered and the two kinds of Rashba spin-orbit 
couplings can drive an ordinary insulating state of the 
silicene to the topological insulator. 



FIG. 4. (color onhne) Phase diagram showing the conduc- 
tance G as a function of A_ri and \r2 for fixing Si — O.lt and 
E^O. 

In order to have a global view on the phase transitions, 
the phase diagram at er^ = O.lt is plotted in Fig. 4. When 
these two kinds of spin-orbit couplings arc tuned prop- 
erly, two kinds of QSHE, QSHEl and QSHE2 appear. 
For QSHEl, the silicene nanoribbon has a large bulk gap 
and there are only gapless edge states in the bulk gap 
because the first spin-orbit coupling can widen the bulk 
gap.^^ On the other hand, for QSHE2, the system has 
a narrow bulk gap and even the gapless edge states is 
located in the bulk band because the second spin-orbit 
coupling can narrow the bulk gap.'^" 
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FIG. 3. The local spin-up (al and bl) and spin-down (a2 
and b2) currents density vector distribution for Si — O.lt, 
Xri = l.Ot and \r2 = 0.2t (QSHEl) with (a) E = O.lt and 
(b) E = -O.lt. 

Next; for QSHEl, the Fermi energy is tuned to E ~ 
O.lt OT E — —O.lt, reaching slightly into the bulk band. 
The configurations of spin-dependent local-currcnt-flow 




FIG. 5. (color onhne) The conductance G as a function of the 
disorder strengths w for (a) Ei — O.lt, \ri = l.Ot and Xr2 = 
0.2t (QSHEl) and (b) Ei = O.lt, Xri = 0.05t and \r2 = 0.6t 
(QSHE2) with the various energy E. The error bars show 
standard deviation of the conductance for 100 samples. 

Finally, we examine the non-magnetic disorder effect 
on the conductance plateau jh of the QSHE. Ran- 
dom on-site potential Wi is added for each site i in the 
central region, where Wi is uniformly distributed in the 
range [— u'/2,w/2] with the disorder strength if. Figs. 5 
(a) and (b) show the conductance G versus the disorder 
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strength at various energy for QSHEl and QSHE2, re- 
spectively. The results show that these quantum plateaus 
are robust against non-magnetic disorder because of the 
topological origin of the edge states. Especially, the 
quantum plateau of QSHEl maintains its quantized value 
very well even when w reaches 2.0t for E = O.Obt, as 
shown in Fig. 5 (a). Wc can also find that the gap- 
less edge states of QSHEl are more insensitive to the 
non-magnetic disorder than those of QSHE2 because the 
bulk gap of QSHEl is larger than that of QSHE2. With 
further increasing of the disorder strength, the conduc- 
tance gradually reduce to zero and the system eventually 
enters the insulating regime. 

In summary, we predict that the QSHE can be induced 
by applying an electric field in silicene even if the intrin- 
sic spin-orbit coupling is very weak. The energy bands, 
the configurations of the spin-dependent local-current- 
flow vector, and the conductance of the system are nu- 
merically studied using the tight-binding Hamiltonian. 
The first and second Rashba spin-orbit couplings, refer- 
ring to the nearest and next-nearest neighbor hoppings 
respectively, can be tuned by the external electric field 
due to the buckled structure of silicene. The staggered 
sublattice potential induced by the external electric field 
open a bulk gap and the gapless edge states are built in 
the gap by the two kinds of Rashba spin-orbit couplings. 
With the help of spin-dependent local-current- vector con- 
figurations, we find that the gapless edge states are in- 
deed spin-filtered. Wc also find that when the two kinds 
of Rashba spin-orbit couplings arc tuned properly, there 
are two types of QSHE, QSHEl with a wide bulk gap 
and QSHE2 with a narrow bulk gap. Moreover, the gap- 
less edge states have also been found to be robust against 
non-magnetic disorder. 
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